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Abstract 

Massive gravity previously constructed as the spin-2 quantum gauge theory is studied 
in the classical limit. The vector graviton field v which does not decouple in the limit 
of vanishing graviton mass gives rise to a modification of general relativity. The modified 
Schwarzschild solution contains a contribution from the t^-field which can be interpreted as 
the dark matter mass density. We calculate the corresponding density profile in the simplest 
spherical symmetric geometry. 
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1 Introduction 



Massive gravity is a controversal subject because it is very difficult to construct it starting from 
a classical lagrangian theory. On the other hand, there exists a powerful method to find the 
quantum theory directly: the derivation of gauge couplings from a cohomological formulation of 
gauge invariance in terms of asymptotic free fields. This method works for massless and massive 
gauge theories equally well. For spin-1 theories this was demonstrated in the monograph [1], 
where the massless spin-2 case is also treated. The massive spin-2 theory was first investigated 
in [2] . The most elegant way to obtain the theory is by assuming the gauge invariance condition 
for all chronological products in the form of the descent equations [3]. These give the total 
interaction Lagrangian including ghost couplings and the necessary coupling to a vector-graviton 
field. By comparison with the couplings derived from expansion of the Einstein Lagrangian in 
powers of the coupling constant we observe that we really have a quantum theory of gravity. 
The massive spin-2 theory corresponds to the Einstein Lagrangian with a negative cosmological 
constant. So this is not directly related to the dark energy (contrary to our earlier belief). 

By analogy with spin-1 one would expect that if the graviton is massive, a gravitational 
Higgs field would be necessary to save gauge invariance in higher orders. This would be a 
natural candidate for dark matter. However, the theory is gauge invariant in second order 
without any additional field [2], and in third order the anomalies checked so far also cancel. A 
gravitational Higgs seems not to exist. 

There is a second option in massive gravity. A massive spin-2 particle has 5 degrees of freedom 
in contrast to 2 of the massless graviton. It turns out that there are even 6 physical modes in 
massive gravity [4]. In order to have a smooth massless limit m ^ of the massive theory, 
one is forced to choose the 6 physical states as follows: two of them are the transversal modes 
of the massless graviton, the remaining four are generated by a vector field (called vector- 
graviton field, see next section) which is required in the massive theory to have it gauge invariant. 
The surprising fact is that for m ^ this vector-graviton field does not decouple from the 
symmetric tensor field /i^*^ which corresponds to the classical g^^'^ of Einstein. Consequently, 
the massless limit of massive gravity is not general relativity; there is a modification due to the 
(now massless) vector graviton field v^. We will show that this modification gives rise to an 
additional force which looks as if it comes from a dark matter density. 

2 Massive quantum gravity 

The basic free asymptotic fields of massive gravity are the symmetric tensor field h^"{x) with 
arbitrary trace, the fermionic ghost u^{x) and anti-ghost u^{x) fields and the vector-graviton 
field v^{x). They all satisfy the Klein-Gordon equation 



(□ + m^)h'^'' = = (□ + m'^)u' 



;^ = (□ -I- m^)u' 



;^ = (□ -I- m^)v^ 



(2.1) 



and are quantized as follows [2] [3] 



[/i-^(x),/i'^^(y)] = -ib'^^t'^D^ix - y) 



(2.2) 



with 



h' 



(2.3) 
(2.4) 




{u^{x),u^{y)} = ir]^^Dm{x-y) 
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[v>'{x),v''{y)]='-'n^^-Dm{x-y), (2.5) 

and zero otherwise. Here, Dm{x) is the Jordan-Pauh distribution with mass m and r^^^^ = 
diag(l, —1, —1, —1) the Minkowski tensor. 

The gauge structure on these fields is defined tlirough a nilpotent gauge charge operator Q 
satisfying 

= 0, = (2.6) 

where O is the Fock vacuum and 

dqh'''' = [Q, h^""] = -^-{d^'u^' + d^v!" - v'"'dau'^) (2.7) 

dQu'' = {Q,u} = 

dQU^' = {Q, u^'] = i{d^h'"' - mv") (2.8) 

dQv'' = [Q,v''] = ~mu''. (2.9) 

The vector-graviton field is necessary for nilpotency of Q. The vacuum CI exists on Minkowski 
background only, therefore, apart from simplicity this background is preferred for physical rea- 
sons. 

The coupling T{x) between these fields follows from the gauge invariance condition [2] 

dQT{x) = idaT"{x) (2.10) 

where T and T" are normally ordered polynomials with ghost number and 1, respectively. In 
addition we may require the descent equations 

^qT" = [Q, T"] = idpT^'f^ (2.11) 

[Q,r"''] =z5^r"^T (2.12) 

where the new T's are antisymmetric in the Lorentz indices. The essentially unique coupling 
derived from (2.10-12) is given by [3] 



-hi^'^Kph + lh^). (2.13) 



Here h = h^'is the trace and a coupling constant is arbitrary. The quartic couplings follow from 
second order gauge invariance and so on. 

We consider the limit m — >^ in the following. The massless limit of massive gravity is 
certainly a possible alternative to general relativity. The new physics comes from the surviving 
coupling term of the vector-graviton 

n = h'^^do^vxd^v^'. (2.14) 
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Since all fields are now massless we may compare with the photon-graviton coupling which 
in linear gravity is given by h°'^tap, where tap is the energy-momentum tensor of the electro- 
magnetic field. It is traceless and conserved, d^tap = 0. Expressing the field tensor F^^'^ by 
the vector potential we get a coupling term of the form (2.14) plus four different additional 
terms. That means Ty is different from photon-graviton coupling, and in fact it does not come 
from a traceless conserved tensor tap- For relativists this may be odd, but we must emphasize 
that our first principle is gauge invariance, general relativity is secondary. 

To understand the situation better we consider a general coupling between h"^ and a sym- 
metric tensor t^'^ with zero gauge variation dqt^'^ = 0: 

T = h'^'^lairta^cVfSu + Va.Vf3i.) + br,ai3V,^.]t'"'. (2.15) 

Prom (2.7) we get 

dgT = -'-{d'^uf^ + - ri'^^dxu^)i2atap + bria^t''^) 

which by partial integration is equal to 

= div + iaul^d'^tap - i{a + b)u^di3t'^ 

where div is a divergence idaT"'- Now gauge invariance (2.10) requires the vanishing of the 
additional terms 

d^ta/3 = (1 + ^)9^tjl. (2.16) 

Using d^v\ = we see that this condition is satisfied by ta/3 = daVxdpv^ (2-14) with a = 2, 
6 = — 1 as it must be, but tag is neither conserved nor traceless. 

To be able to do non-perturbative calculations in massive gravity we look for the classical 
theory corresponding to the couphng (2.13). It was shown in [1] sect. 5.5 that the pure graviton 
couplings hdhdh correspond to the Einstein-Hilbert Lagrangian 

i^EH = ^V^R (2.17) 
in the following sense. We write the metric tensor as 

^g^''^ = ^^'^ + ^ih^'^ (2.18) 

and expand Leh in powers of k. Then the quadratic terms 0{k^) give the free theory, the cubic 
terms 0{k^) agree with the pure graviton coupling terms in (2.13) up to a factor 4 and so on. 
The mass dependent terms are obtained if one introduces a cosmological constant —A = m^/2 
in (2.17) [3], the minus sign is due to the different convention in astrophysics [8]. As already 
said we consider the massless limit in the following. 
To obtain Ty (2.14) we must add 

Lv = ^g'^^davxdpv^ (2.19) 

to Leh- One may be tempted to write covariant derivatives V a instead of partial derivatives 
in order to get a true scalar under general coordinate transformations. But this would produce 
quartic couplings containing v\ and such terms are absent in the quantum theory ([2], eq.(4.12)). 
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For the same reason the Lorentz index A in v\ is raised and lowered with the Minkowski tensor 
rffj,,^, but all other indices with ^^j,. Again this might look odd to relativists, but the correct 
classical expression is dictated by the quantum theory not by something else. Note that the 
gauge variation (2.7) is just the infinitesimal form of general coordinate transformations, so that 
this symmetry is built in on the quantum level not on the classical one. 
The total classical Lagrangian we have to study is now given by 

Ltot = Lg + Lm + Ly (2.20) 

where 

we have introduced the velocity of light and Newton's constant for later discussions. Lm is given 
by the energy-momentum tensor of ordinary matter as usual. Strictly speaking we should 
also add a term for the ghost coupling. At present we leave this out because probably these 
unphysical degrees of freedom don't contribute in the classical theory. 

Calculating first the variational derivative with respect to g^i, we get the following modified 
Einstein equations 

1 „ SttG^ 167rG/ ^ A 

Rij.u - -^9iivR = —^\-diJ.v\OvV^+ 

+\g^ug''^davxdpv'^). (2.22) 

Secondly the variation with respect to v\ gives 

daiV^g'^^ds.v^) = 0. (2.23) 

After multiplication with \j ^—g this is the Laplace-Beltrami or rather the wave equation in 
the metric 5"^. The coupled system (2.22) (2.23) is the modification of general relativity which 
we want to study. 



3 Modified Schwarzschild solution 

Wc wish to construct a static spherical symmetric solution of the modified field equations. 
Following the convention of [5] we write the metric as 

where v and A are functions of r only. We take the coordinates = ct, = r, x"^ = x^ = ip 
such that 

500 = e", gu = -e^ 

522 = -r^, gss = -r"^ s'm^ ^ (3.2) 

and zero otherwise. The components with upper indices are the inverse of this. The determinant 
comes out to be 

g = det£?^^ = -e'^+^r^ sin^ (3.3) 
Next we compute the Christoffel symbols 

J^oo — 2 ' ^00—2 ' ^22— 
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We obtain 



rl^ = -re-^sm^^, Tj^ = K rig = - sini?cosi?, rfg = ^, rig = cot 7? 

and zero otherwise, the prime denotes the derivative with respect to r always. Prom this the 
contracted curvature tensor is given by 

f)-pX /5pA 
R _ fMi' _ " t^X I pA _ -pQ pA 

'''' ~ dx'' 

B^^ = \e'-\-" + \-"-\-'>^ + -/) (3.4) 

Ri, = -\u" + \\'{v' + -^)-\v'^ (3.5) 

R22 = e-\'-X' -'-u'-\) + l (3.6) 

R^z = sin^ i»R-22. (3.7) 

For the scalar curvature we then find 

R = e-Hv" + 1^.'^ - \v'>! + 2- -2^ + ^)-^. (3.8) 
Now we are ready to calculate 



G; = i?;-^<^;:i?, (3.9) 

G'o=e-\^-^)A (3.10) 



It 12 I \i \i I 



r 4 4 2r 2r 
The energy momentum tensor of ordinary matter is assumed in the diagonal form 

= diag(Qc2, -pi, -p2, -Ps) (3.13) 

where q{r) is the (inertial) mass density and Pj{r) the presure components. 

After collecting these well known results [5] we turn to the new contributions involving the 
vector-graviton field v'^. Since 

g(i/+A)/2 q2 



sin ■& dip^ ' 

the Laplace-Beltrami equation (2.23) can be written in the form 



(3.14) 



f (e('^-^)/V2f)z;« = e(^+^)/2L\„, (3.15) 
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where is the quantum mechanical angular momentum operator squared. Consequently the 
angular dependence of is given by spherical harmonics Yj"^('!?, </?): 

v,{r,i^,^)=v^{r)Yri^,'p) 
and the radial function ^^(r) satisfies the following radial equation 

|.(e(^-A)/V^^)= + l)e(^+^)/^..(r). (3.16) 

For what follows we consider the case I = where we get the simple result 

<(r) = ^e(^-)/^ (3.17) 

here C^, k = 0, 1, 2, 3 are constants of integration. 

Now wc are ready to write down the modified radial Einstein equations. In (2.22) we raise 
the second index 

G'/T = + V'^d-vxdpv^ - 2d^v,d''v'^). (3.18) 

G° = e-^(^ - = ^(g(r)c + ^o(r)) (3.19) 



. . V - ^ - ^ . ^^(^ifi . (3.22, 

where wq,wi,W2,W3 are the contributions from the f -field in (3.18). The two equations (3.19) 
(3.20) together with (3.15) are sufficient to calculate all unknown functions A, v, v^. That means 
the equations (3.21) (3.22) are additional constraints. It is a delicate problem to satisfy them. If 
we assume p2 = Ps, we must have W2 = ws and this is only possible for I = 0. Then v^, depends 
on r only. We restrict ourselves to this case in the following, i.e. pj = p. 
From (3.17) we now obtain the following contributions of the ?;-field 

g^'div^div- = ^e-^(Cf + C| + C| - C^) 

= ^e-'' (3.23) 
(J 

Wo = —7e~'' = W2 = W3 (3.24) 
Wi = —wq. 

We see in (3.19) that the mass density q{r) gets enlarged by the quantity WQ{r)/c which we 
shall call dark density. We omit the notion "matter" because this energy density comes from 
the vector-graviton field, it is a relic of the massive graviton. 



and zero otherwise, so that 
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Next by suitable combination we simplify the equations. Adding (3.19) to (3.20) we get 

^(A' + ^^')=5(9 + ^) (3.25) 

where 

SttG 

9 = 3.26 
Eliminating p from (3.20) and (3.21) we obtain 

u" = 2e\-- + 2g^)+—- — + + ^. 3.27 

Next we differentiate (3.20) with respect to r and use (3.27): 

- -e-^f (V + u') + 4^wo. (3.28) 

Vc^ c ^ 2r cr 

Substituting (3.25) inhere we finally arrive at 

^ = -^{1 + ^ + 1^0) (3.29) 

where (3.24) has been used. This differential equation for the pressure will be used instead of 
the second order equation (3.21) in the following. 

For a first orientation we want to solve the equations neglecting the ordinarary matter, i.e. 
q = = p. Then we find from (3.25) A' + i/' = or 

u{r) = Di- A(r) (3.30) 

where Di is a constant of integration. Using this is (3.24) we have 



wo = §e\ (3.31) 



Multiplying (3.19) by 

we write the equation in the form 

where 



y = re~^ (3.32) 



y' = l- Stt— (3.33) 
ry 

D, = (3.34) 



Although this equation (3.33) looks rather simple, we did not succeed in expressing the solution 
by known special funtions, but a numerical solution can be easily obtained. The dark density 
wq (3.31) is given in terms of y (3.32) by 

^0 = ^-TT- (3-35) 
For large r we find the following power series expansion 

y = r-ai + — + —^ + ... (3.36) 
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where 

a2 = 87rr>3, 03 = 47rD3ai,--- (3.37) 
ai is related to the total mass M by the Schwarzschild relation 

ai = ^ . (3.38) 

This can be taken as the start of the numerical integration of (3.33) from large to small r. Note 
that the Schwarzschild horizon disappears, but there remains a singularity at r = 0. It is of the 
form 

y = 4, -ttDs logr + a + -r^—- 
y logr 

This is too weak to render v^wq (3.35) integrable at r = 0. We must conclude that the classical 
theory breaks down at r = 0. This is a phenomenon known from ordinary Einstein equations. 

Without better options we cut off the singularity by assuming wq = const, inside a core 
radius Vc, for r > rc wq is given by (3.35). The core radius Vc is determined as follows: We 
compute the total mass 



M = —( f + I )r^u;o(r) 



CO 







CO 



c V 3 J ryJ 



= -f4^ + ^ l{i-y')dr 

c V3y(re) BttDs J 
An D2 (? 

where (3.34) has been used. By (3.38) M drops out and we get Tc as the zero of the equation 

rc-y{Tc) = ^-^^y (3.39) 

3 y{rc) 

Varying this enables us to vary Tc and the radius of the halo, if we keep its total mass M, 
i.e. a\ (3.38) fixed. The differential equation (3.33) is invariant under the scale transformation 

r ar, y ay, ci^D^, ai aai. (3.40) 

This implies that the length scale of the halo can be freely adjusted. 

In astrophysical applications we measure r in kpc. Let us consider a halo of 10^"^ solar masses, 
then ai (3.38) is equal to 10""^ (in kpc). Taking D3 of the order 10~^ we find rather small core 
radii of a few kpc. Precisely, our numerical results agree with the linear relation 

rc = 3352^)3. (3.41) 

Outside rc, y is given by the four terms in (3.36) up to promille accuracy. That means the dark 
density profile (3.35) can be approximated by 

^0 = — ^ 2T ^ ^ (3.42) 

r — air"^ + a2r + 03 



00 
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Then the dark mass interior to r is given by 

r r 

wor'^dr = ATrD2 / 5 . (3.43) 

J - air^ + a2r + 03 

Expanding this for small r (but > Tc) we find a linear dependence 

M (r) = Qo + air + . . . (3.44) 

This corresponds to constant orbital velocity as observed in dark halos of spiral galaxies. 
On the other hand Navarro, Frenk and White [6] have found the following profile 

g{r) 



from N-body simulations, where is the so-called scaling radius and Qs the scaling density. 
In contrast to our result (3.42), g{r) only decreases as for large r. But then r'^Q{r) is 
not integrable at infinity. A profile with tail was proposed by Hernquist [7] for elliptical 
galaxies. It remains to be seen which profile gives the best agreement with reality. 
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